Passive Scalar Evolution in Peripheral Region 
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We consider evolution of a passive scalar (concentration of pollutants or temperature) in a chaotic 
(turbulent) flow. A universal asymptotic behavior of the passive scalar decay (homogenization) 
related to peripheral regions (near walls) is established. The passive scalar moments and its pair 
correlation function in the peripheral region are analyzed. A special case investigated in our paper 
is the passive scalar decay along a pipe. 
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INTRODUCTION 

The problem of mixing attracts a great attention both 
due to its fundamental significance and a variety of ap- 
plications. Mixing rate of an additive to a fluid is very 
sensitive to the character of the hydrodynamic flow ex- 
cited in the fluid. Recently an essential breakthrough 
was achieved in the theory of the so-called passive scalar 
(which can be concentration of pollutants or tempera- 
ture) in chaotic and turbulent flows (see reviews Q>Q|)- 
The word "passive" means that a feedback of the scalar 
field to the flow is negligible. It is correct for dilute so- 
lutions of pollutants and for relatively weak fluctuations 
of temperature. The passive scalar in a fluid is subjected 
to diffusion (thermodiffusion) and advection, leading to 
an evolution of its spatial distribution. A final result 
of the evolution is a homogeneous state of the passive 
scalar. We are interested in the passive scalar decay, 
that is in the laws governing the homogenization of the 
passive scalar field in chaotic and turbulent flows. 

As it is noted in the paper at large Peclet number 
Pe the homogenization of the passive scalar field in the 
peripheral regions is slower than in the bulk since mixing 
is suppressed near walls of the vessel. Then the advanced 
stages of the passive scalar decay are dominated just by 
the peripheral regions. In the paper Q main stages re- 
lated to the passive scalar evolution in the peripheral re- 
gion were established (see also the paper 0). Here we de- 
velop further the theory of the passive scalar evolution in 
the peripheral region having in mind two different phys- 
ical situations. The first case is the passive scalar in the 
viscous boundary layer of the developed (high-Reynolds) 
turbulence (explanations of its properties can be found in 
the book of Monin and Yaglom j5j). The second case is 
a chaotic flow. A perfect example of such flow is the so- 
called clastic turbulence state revealed by Groisman and 
Steinberg in dilute polymer solutions. Both the veloc- 
ity in the viscous boundary layer (in the first case) and 
the chaotic velocity (in the second case) can be treated 
as smooth. The smoothness of the random velocity field 
is one of key ingredients of our analysis. 

The theory of the passive scalar in a smooth random 
flow was pioneered by Batchelor and Kraichnan Q 
who considered the cases of velocities long-correlated and 



short-correlated in time, respectively. The case of arbi- 
traryvelocity correlation time was considered in the pa- 
per In the works the passive scalar was assumed to 
possess a stationary statistics due to pumping. Prop- 
erties of the passive scalar decay in an unbounded sta- 
tistically homogeneous turbulent media are also known. 
As it was demonstrated in the paper |1C| . in the iner- 
tial range of the developed turbulence the scalar decays 
in accordance with a power law. Besides, it is interest- 
ing to consider the passive scalar decay in the Batchelor 
region of scales (below the Kolmogorov viscous scale), 
where the velocity field is smooth (for details see, e.g., 
the book of Batchelor 11]). The decay law of the passive 
scalar in the region is exponential, as it is demonstrated 
in the papers [2> El- The same consideration is ap- 
plicable to the elastic turbulence. In the paper Q the 
combined case was considered, when both the inertial re- 
gion of scales and the Batchelor region of scales are taken 
into account. Then the scalar decay is dominated by ed- 
dies from the inertial interval (excluding, maybe, some 
initial stage of evolution) and is, consequently, governed 
by a power law. 

As it is noted in the paper [j| the theory of the pas- 
sive scalar in the random flow, developed for the bulk, 
needs a modification for the peripheral region even for 
the smooth flow. The reason is that the approximation 
of the velocity by linear profiles, used for the bulk, fails 
in the peripheral region. However, near a wall, the ve- 
locity possesses a definite dependence on the separation 
from the wall, explaining the universal properties of the 
passive scalar decay in the peripheral region. Another 
property, simplifying an analysis, is the slowness of the 
passive scalar mixing in the peripheral region, enabling 
one to treat the velocity as short-correlated. The first 
problem, we consider, is an evolution of the concentra- 
tion of pollutants, when the boundaries are considered as 
impenetrable for the pollutants. This approach can be 
extended to the case of the binary chemical reactions (see 
the paper £|). Another problem, which can be treated 
inside our approach, is the temperature relaxation in the 
bulk, if the temperature at the walls of the vessel is fixed. 
Then a heat flow from the boundary to the bulk is forced 
which that is governed just by the velocity fluctuations 
in the peripheral region. The passive scalar evolution can 
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be examined also in a fluid (where random motion is ex- 
cited) flowing through a pipe. Then the role of time is 
played by the coordinate along the pipe. This setup is 
closely related to the experiment by Groisman and Stein- 
berg 0. 

A remarkable property of a turbulent flow is its strong 
intermittency leading to the so-called anomalous scaling 
of the velocity correlation functions, which are dependent 
on the integral scale of turbulence at a power, referred 
to as the anomalous scaling index (see, e.g., the book of 
Frisch (l5||). Theoretically, an existence of the anoma- 
lous scaling of such kind was established for the pas- 
sive scalar in the framework of the so-called Kraichnan 
model (see the papers ll^j). One could anticipate that 
the same phenomenon takes place for the passive scalar 
decay. We demonstrate that in the peripheral region the 
passive scalar moments possess, indeed, an anomalous 
scaling, which is in some sense extreme. Namely, all the 
moments of the passive scalar damp in accordance with 
the same power law at increasing the separation from the 
wall. 

The paper is organized as follows. In Section [I] we 
present general equations, which describe the passive 
scalar evolution in the peripheral region. In Section [H] 
we analyze the passive scalar evolution near the wall of 
the vessel. In Section UTTl we consider the passive scalar 
decay along a pipe. Some general remarks and a short 
comparison with experiment are presented in Conclusion. 

I. GENERAL RELATIONS 

Advection of a passive scalar field 9 by a moving fluid 
(accompanied by the passive scalar diffusion) is described 
by the equation 

d t e + vve = nv 2 e , (1.1) 

where v is the flow velocity and k is the diffusion coef- 
ficient. Below, the fluid is assumed to be incompressible 
(that is Vt> = 0). A formal solution of the Cauchy prob- 
lem for the equation can be written as 

9(h) = Tex p{^ 2 dt + kV 2 ]| 9(h) , (1.2) 

where Texp means a chronologically ordered exponent. 
Of course, some boundary conditions for the passive 
scalar 9 should be introduced. There could be two differ- 
ent types of the boundary conditions. If 9 is temperature 
(and walls are made of a well heat conducting material) 
then 9 is fixed at the boundary. If 9 is density of pollu- 
tants and the boundary is impenetrable for the pollutants 
then the gradient of 9 in the direction perpendicular to 
the boundary is zero (that corresponds to zero pollutant 
flux to the boundary). 

We consider a random flow which has to be character- 
ized statistically: via correlation functions. The correla- 
tion functions are averages over time, they can be treated 



also as averages over velocity realizations. The flow is as- 
sumed to be statistically homogeneous in time, whereas 
there is no homogeneity in space (because of the bound- 
ary effects). Here we consider the flow in a closed vessel 
where the average velocity is equal to zero. A generaliza- 
tion to the case when the average velocity is non-zero is 
trivial (it can be found in Section lTTTl where a fluid pushed 
through a pipe is treated). The pair velocity correlation 
function (v a (h, ri)vp(t 2 , r 2 )) depends on the time differ- 
ence h — h only (due to the assumed time homogeneity), 
and on coordinates of both points r\ and r 2 . 

We have in mind two different situations. First, it 
is the viscous boundary layer of the developed high- 
Reynolds turbulence (for details see, e.g., the book 0). 
Second, it is the peripheral region of a chaotic flow. As a 
perfect example of the chaotic flow, the elastic turbulence 
can be noted [6j . In both cases the velocity is smooth in 
the peripheral region we are interested in. However, the 
velocity field possesses some peculiarities, related to zero 
value of the velocity at the boundary. That makes the 
passive scalar decay in the peripheral region slow. Since 
the velocity correlation time is determined by dynamics 
in the bulk, which is relatively fast, at examining the 
passive scalar evolution in the peripheral region the ve- 
locity can be treated as short correlated in time. It is well 
known, that in this case closed equations for the passive 
scalar correlation functions can be derived. Below, we 
demonstrate principal steps of the derivation, based on 
Eq. (Ol) . 

Let us examine the passive scalar evolution on a time 
interval (h,t 2 ) taking the difference h — h much larger, 
than the velocity correlation time r, but much smaller, 
than the characteristic mixing time (the gap between the 
mixing time and the velocity correlation time exists due 
to the noted weakness of mixing in the peripheral region) . 
The last condition enables one to produce an expansion 
of the T-exponent in Eq. (|1.2|) . One may keep only two 
first terms of the expansion: 

9(h) « 9(h) + (h - h)nV 2 9(h) - f 2 dtv(t)V9(h) 

Jti 

+ f 2 dt f dt'v(t)V [v(t')V9(h)} ■ (1.3) 

The next step is in averaging over the velocity statistics 
inside the interval (ti, t 2 )- This averaging is independent 
on the velocity profiles at t < h and t > h due to the 
condition h — h 3* T - Averaging the expression (|1.3fl . 
one obtains a relation for the average value (9) of the 
passive scalar 

(9(h,r)) - (9(h,r)) = (h - h)^ 2 (9(h, r)) 

+ {h - h)V a [D aP (r,r)V p(0(h,r))] , (1.4) 

D ad {rx,r 2 )= / dt(v a (t, ri )v 3 (0,r 2 )) (1.5) 
Jo 

where a fast enough decay of the pair velocity correla- 
tion function with t is implied. We used also the incom- 
pressibility condition Vv = 0. The quantity (9) means a 
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value, averaged over the velocity fluctuations. The quan- 
tity D a p(r,r), entering Eq. I|1.4|) . is usually called the 
turbulent diffusion tensor. If t 2 — ii is smaller than the 
mixing time, then the right-hand side of the equation 
(|1.4() is a small correction to (9) . Therefore the equation 
(|1.4|) can be rewritten in the differential form: 

dt(0) = V Q [D a p{r, r)V p (9)} + nV 2 {9) . (1.6) 

Analogously, starting from Eq. (|1.2I) . one can derive 
closed equations for high-order correlation functions of 9. 
Say, the equation for the pair correlation function F is 

8 t F(t,r 1 ,T 2 )=K(yl + Vl)F (1.7) 
+V iQ [^(r 1 ,ri)V 1 /]+V 2a [D a0 (r 2 ,r 2 )V 20 F] 
+V iQ {D a0 (r 1 ,r 2 )V 2l3 F}+V 2a [D a0 (r 2 ,n)Vi F] , 
F(t,r 1 ,r 2 ) = (9(t,r 1 )9(t,r 2 )). (1.8) 

Generally, the equation for the n-th order correlation 
function F n of the passive scalar is 

n 

dtF n = K ^ VliFn 
m— 1 

n 

+ Vma [D a p(r m ,r k )V k0 F n ] , (1.9) 

m,k—l 

F n (t, ri ,...,r n ) = (6{t, ri )...9(t,r n )). (1.10) 

The structure of the equation i|1.9fl is transparent: the 
evolution of the passive scalar correlation function is de- 
termined by the molecular diffusion (the first term in the 
right-hand side) and by the turbulent diffusion (the sec- 
ond term in the right-hand side). 

Note that if the diffusivity is negligible then it is 
possible to obtain a closed equation for the moments 
(9 n {r)) = F n {r, . . . , r) from Eq. lfO?|l 

d t (9 n (r)) = V Q [D a0 (r, r)V (O n (r))] , (1.11) 

which is identical to the equation (|1.6(l . without the 
molecular diffusion term. Eq. I|l.ll(l is a direct conse- 
quence of the relation dt{9 n ) = —vW(9 n ), which follows 
from Eq. Ijl.lfl . if the molecular diffusion is neglected. 

The turbulent diffusion tensor D can be estimated as 
D ~ V 2 t, where V is characteristic value of the velocity 
fluctuations, and r is the velocity correlation time, that 
can be estimated as r ~ L/Vl- Here L is the size of the 
viscous boundary layer in the case of the high-Reynolds 
turbulence and the size of the vessel in the case of the 
chaotic flow (elastic turbulence), and Vl is characteristic 
value of the velocity fluctuations in the bulk. Therefore, 
the weakness of the passive scalar decay is explained by 
smallness of the ratio V/Vl inside the peripheral region. 
Note that for the elastic turbulence r is determined by 
the polymer relaxation time and the condition r ~ L/Vl 
is no other than the Lumley criterion of strong polymer 



elongation formulated in the papers [iTj (see also the 
papers (lSjl. 

Let us explain physical meaning of the passive scalar 
correlation functions. They are averages over the velocity 
statistics. Therefore, to obtain the correlation functions 
experimentally or in numerics, one has to measure the 
passive scalar decay many times (for many realizations 
of the velocity field) and then to average the result over 
the attempts. Initial conditions for the passive scalar 
field are implied to be fixed at the averaging procedure. 
If we consider the case of a fluid pushed through a pipe, 
then the passive scalar correlation functions are station- 
ary. Then they can be treated as averages over long time. 

Below we consider a passive scalar evolution in the 
peripheral region. We assume that mixing already pro- 
duced the homogeneous distribution of the passive scalar 
in the bulk (remind that the process in the bulk is much 
faster than in periphery). And we subtract from 9 a con- 
stant, corresponding to the bulk value of 9 (this redefini- 
tion does not change the equations describing the passive 
scalar evolution). By other words, the value of 9 tends 
to zero when we go away from the boundary. Next, the 
homogenization of the passive scalar along the boundary 
is much more effective than its homogenization in the di- 
rection perpendicular to the wall. Therefore, analyzing 
the evolution of the passive scalar after the homogeniza- 
tion in the bulk, we can treat its initial condition 9q as 
dependent mainly on the separation from the wall and 
practically independent of coordinates along the wall. 



II. DECAY IN VESSEL 

Here we consider the passive scalar decay in a closed 
vessel. If the walls of the vessel are smooth and their cur- 
vature is of the order of the vessel size, then at considering 
the peripheral region the boundary can be treated as flat 
in the main approximation. Then it is possible to intro- 
duce the orthogonal reference system, where the coordi- 
nate q measures separation from the boundary and are 
coordinates along the wall, the subscript /j (v) running 
over 1 and 2. In the reference frame the incompressibil- 
ity condition is written as d q v q + dvi/dri + dv 2 /dr 2 = 0, 
where v q is the velocity component perpendicular to the 
wall and are the components along the wall. Since 
tends to zero as q — > then oc q near the bound- 
ary, and incompressibility leads to the proportionality 
law v q oc q 2 . This is the main feature of the velocity in 
the peripheral region. 

For the flat wall, it is naturally to assume homogeneity 
of the velocity correlation functions along the wall and 
also their isotropy in the planes parallel to the wall. Then 
we obtain a general expression for the components of the 
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turbulent diffusivity tensor l|1.5fl 



A. Average Scalar 



D^{r u r 2 ) = H 1 (g)5^ l/ q 1 q 2 + H 2 (g)g^g u q 1 q 2 , 
D qu {r 1 ,r 2 ) = -\H[( e )^qlq2 



"2 [e H L(Q) + SHzie)] Qvq\qi-, 



V m {r u r 2 ) = \lI> x {Q) e -±q xq l 
z g 



1 



[qH' 2 {q) + 3H 2 (g)} g^q 2 , 



2 

D qq (r u r 2 
1 



me) 



h'({q) 



2 2 
1ll2 



^ [6H 2 (g) + 6gH 2 (g) + g 2 H! 2 '(g)] qfq 2 , (2.1) 

where g^ = rip — r 2{i and H±, H 2 are some func- 
tions of g. The expressions (|2.1f) satisfy the conditions 
Vi a D Qj g(ri,r 2 ) = = ^7 2 pD af3 (r 1 ,r 2 ), following from 
the incompressibility condition. 

Since we consider the peripheral region where the ve- 
locity field is smooth, then both Hi and H 2 have regular 
expansion in g (containing even powers) at small g, that 
is at g <C L (remind that L is the thickness of the periph- 
eral region which is the thickness of the viscous boundary 
layer in the case of the developed turbulence and is of 
the order of the vessel size for the elastic turbulence). In 
this limit the main contribution to the turbulent diffusion 
tensor is determined by first terms of the expansion: 



H w -(» + 3H 20 /2)g 2 
Duu ~ HioSn^qiq 2 , D qq 



H 2 

mill- 



fla>,(2.2) 
(2.3) 



The quantity fj, characterizes the flow intensity near the 
boundary. It can be estimated as fi ~ Vl/L 3 , where, as 
previously, Vl is the characteristic velocity fluctuation 
in the bulk. In the framework of the Karman-Prandtl 
theory of the viscous boundary layer (for details see, e.g., 
the book p|) the width of the layer is estimated as L ~ 
v/Vl (where v is kinematic viscosity of the fluid) and 
one finds /i ~ V^v~ 3 . For the elastic turbulence L is the 
vessel size and /i ~ V^^~ 3 Re~ 3 , where Re = VlL/v is 
the Reynolds number. 

Comparing the turbulent diffusion term for the mo- 
tion perpendicular to the wall, which can be estimated 
as /j,q 4 d 2 , and the molecular diffusion term nd 2 , one finds 
the width of the boundary diffusion layer 



[k/h) 



1/4 



(2.4) 



We assume that ru <C L. The relation can be rewritten 
as Pe 3> 1 where Pe is the Peclet number: Pe = VlL/k. 
As it follows from Eq. L/r bl ~ Pc 1/4 . For the 

case of the viscous boundary layer, one obtains Pe ~ Sc, 
where Sc = vjn is the Schmidt number. For the elastic 
turbulence, one finds a slightly different estimate Pe ~ 
Sc • Re. 



Here we consider the simplest possible correlation func- 
tion: the average value of the passive scalar (8). The 
quantity essentially depends on the separation from the 
wall q and slowly depends on the coordinates along the 
wall. We ignore the last dependence, assuming that its 
characteristic length is the vessel size. Then we find from 
Eqs. I|1.6I2.3|) the following equation for the object 



d t (9) = ^d q ( q i d q (e)) + Kd 2 (e). 



(2.5) 



The advection term (with /x) in Eq. I|2.5|) dominates at 
q 3> ?'m and the diffusion term (with k) dominates at 
g C fa, where ru is thickness of the boundary diffusion 
layer defined by Eq. I|2.4I) . As we explained in the last 
paragraph of SectionQJ at solving the problem we should 
assume lim^oo (9) — 0, since large q correspond to the 
bulk, where 9 is assumed to tend to zero. The condition 
is implied below. 

We examine the passive scalar evolution in the periph- 
eral region, which begins after that its homogenization 
in the bulk is finished. Then the initial distribution of 
the passive scalar has the characteristic length L. The 
subsequent evolution is divided into two stages. At the 
first stage the thickness S of the layer, where 9 is con- 
centrated, diminishes as S — (/it) -1 / 2 . When S reaches 
ru, the second stage starts, which is characterized by the 
fixed spatial scale ru- 

At treating the first stage one can omit the diffusion 
term (with k) in Eq. I|2.5|l . that leads to the equation 



d t {9) =^ q (q%(9}) 



(2.6) 



Looking for a solution (9) = exp(— st)tp s (q), one obtains 
from Eq. 



7T / S 

2 W. 



Using the orthogonality relation 

7TS 3 / 2 



v 3/4 

j J3/2 (V^M'T 1 ) • (2-7) 



dqtp s (q)tp a (q) 



,1/2 



10 M 
one finds a general solution of Eq. I|2.6|l 

st r°° 



5(s-a) , 



(0(t,q)) = / dstfisiq) 



y/JIe- 



,3/2 



dq' ip a (q')0 o (q') (2.8) 



in terms of the initial passive scalar distribution. The 
approximation H2.8(l is correct provided S ^> ru (then it 
is possible to neglect the diffusion boundary layer where 
diffusion is relevant). 

We start from a distribution of 9 with the characteristic 
length L. Therefore at times, when 5 <C L and at q -C L 
one can substitute 9(t = 0) in Eq. (|2.8|) by fl = 9(t = 
0,q = 0). Taking then the integrals in Eq. (|2.8|) . one 
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derives 



Wt,g)> 



#0 



ds 



'7T (J 



■ exp 



(2.9) 



This profile has a universal form insensitive to details of 
the initial distribution of 9. If q S then one obtains 
fromEq. $2~ty 



(2.10) 



If q <C S then we find (9) = i?o . So, the value of 9 is 
practically unchanged inside the layer q < S. 

Note that though the equation ((2.61) has form of the 
conservation law for (9), the total amount of the pas- 
sive scalar in the peripheral region J dq (9) appears to 
be time dependent, if the expression l|2.9|l is integrated. 
The reason is that the considered solution corresponds 
to non-zero passive scalar flux directed to large q, i.e. to 
the bulk, which can be treated as a big reservoir. This 
flux nq A dq(9) can be obtained directly from the asymp- 
totic expression l(2.1(J() . Note also that the passive scalar 
evolution at the first stage is insensitive to the boundary 
conditions, and therefore it is described identically for 
the concentration of pollutants and temperature. 

Now we analyze the passive scalar behavior at the sec- 
ond stage, then the diffusion term can not be ignored. Let 
us first examine the case when the passive scalar repre- 
sents the concentration of pollutants. Then the equation 
(|2.5|) has to be supplemented by the boundary condition 
d q {9) = at q = 0, which means zero flux of pollutants 
to the boundary. At long times, only the contribution 
related to the minimal (by its absolute value) eigen value 
of the operator in the right-hand side of Eq. (|2.5|l is left. 
That leads to the exponential decay (9) cx exp(— jt), the 
decrement 7 is related to the energy of the ground state. 
The value of 7 can be found numerically, it is 7 = Cev/k/I, 
ce ~ 1.81. The asymptotic behavior of (9) can be related 
to the initial value of the passive scalar $0 near the wall: 
(9(q = 0)) = c $o ex P( — it), where c ~ 1.55. The to- 
tal amount of the scalar near the boundary behaves as 
J dq (9) = cid^ru exp(— jt), where c\ « 1.55. 

Let us now consider the case when the passive scalar 9 
represents temperature, assuming that it is fixed at the 
boundary: 9(q = 0) = $q. Then after the first stage a 
quasi-stationary distribution of (9) is formed, since the 
bulk can be treated as a big reservoir having a constant 
temperature. This quasi-stationary distribution can be 
found directly from the equation l|2.5|) where the term 
with the time derivative has to be omitted: 



2V2 



dqi 



[iq\ + K 



(2.11) 



At q 3> ru we find, again, (9) oc q~ 3 . That corresponds 
to a non-zero passive scalar flux (heat flux) to the bulk. 
This flux is time-independent in the case considered. 



B. High moments 

As we already noted, at the first stage the diffusive 
term in the equation for the passive scalar correlation 
functions can be omitted. Then the closed equation 
(|1.11|) for high moments of the passive scalar is cor- 
rect. The high moments (similar to the first one) de- 
pend mainly on q and we ignore their dependence on the 
coordinates along the wall. Then, substituting into the 
equation the expressions 1(2.3(1 . one finds 



d t {9 n )= f xd q (q%{9 n )) 



(2.12) 



which is a generalization of Eq. 1)2.6(1 . Its solution can 
be written as Eq. (|2.8|) 



dsip s (q) 



,3/2 



dq' ip s (q')9Z(q'). (2.13) 



The expression (12.13(1 implies that 5 3> ru since only at 
this condition it is possible to neglect the region q ~ ru 
where the diffusion is relevant. Since the initial distribu- 
tion of 9 has the characteristic length L, at the condition 
S<Iwc obtain, again, a universal expression 



'■K q 



■ exp 



"V 



(2.14) 



which is a generalization of Eq. 1(2.91) . The expression 
shows that in the region q > 5 (9 n ) m tf%5 3 /(6y/nq 3 ). 

Actually, the expression 1(2.13(1 is correct for any aver- 
aged local function of the scalar 9. We can use it to find 
the local scalar PDF P(t, q, 9) = (6 [9 - 9{t, q)]). Let us 
first rewrite the expression like ((2.13(1 as 

1 r°° r°° qV 

xJ 3/2 (k/q)J 3/2 (k/q')P(0,q'), (2.15) 
One can take the integral over k in Eq. 1(2.15(1 to obtain 




7 exp 



Aq 2 nt Aq l2 nt 



|cosh 


1 




2qq'fit 



2[itqq sinh 



2qq' fit 



P(0,q') 
. (2.16) 



Let us assume that the initial scalar distribution is a 
monotonic function 9o(q), which is equal to zero at q — > 
00 and reaches a maximum value at q — > 0. Then P(t = 
0,q,9) = -[l/9' (q)]S[q-q (9)], where the function q (9) 
is determined from the relation #0(90) — 9. Calculating 
the integral 1(2.16(1 one finds 



P(t,q,9) = 



1 



qqo\9' (q )\ 



g [ - - — ] (1 - 2[itqq ) 

q go/ 



-9 ( - + — ) (1 + ^tqq ) 

\q qo 



(2.17) 



1 



: exp 



X 
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When t grows, the boundary layer, determined by 6, 
shrinks, that is the characteristic value of q decreases. 
Besides, qo is fixed at a given 0. Therefore one can ex- 
pand the expression in the square brackets in Eq. (|2.17|) 
in 1/qo to obtain a universal probability distribution 

p = iWPFW^ exp (" v^) ' (2 ' 18) 

Unfortunately, the expression (|2.18l) cannot be utilized 
for calculation of the moments of 0, since the integrals 
J d0 n P(0) diverge near the maximum value of 0, i?o- 
The reason is that the expression l|2.18|l is correct only if 
qo(0) 3> which is violated at small qo, corresponding 
to close to $o- 

At the second stage diffusion starts to be relevant, and 
it is impossible to obtain closed equations for the mo- 
ments (0 n (t,q)}. To find the moments one has to solve 
the complete equations <|1.9|l for the passive scalar cor- 
relation functions, that is a complicated problem. One 
can say only that due to linearity of the problem, the 
asymptotic in time behavior of the correlation functions 
(and, consequently, moments) is determined by the min- 
imal eigen value of the operator in the right-hand side 
of the equation l|1.9l) . Therefore for the concentration of 
pollutants (0 n (t,q)) cx exp(— j n t) where 7„ ~ y/Kfi. The 
estimate can be obtained by equating the time derivative, 
the turbulent diffusion term and the molecular diffusion 
term in Eq. (|1.9() . If we consider the region q ru then 
it is possible to neglect the molecular diffusion term in 
Eq. H1.9|l and we return to the closed equation l|2.12|l . 
There the time derivative d t can be substituted by — j n 
and we conclude that at the same condition q 3> ru it is 
possible to neglect the term with the time derivative in 
Eq. (|2.12|l . Therefore we obtain, again, the asymptotic 
behavior 

(e"(f,?))«f 3 , (2.19) 

at q 3> ru- The situation here is analogous to one for 
the average passive scalar, since the equation (|2.12|) has 
the form of conservation law for the high moments. The 
laws (0 n (t,q)) cx q~ 3 correspond to non-zero fluxes of 
high degrees of the passive scalar. 

Now we can turn to the case when the passive scalar 
is temperature, fixed at the boundary. If the temper- 
ature in the bulk is different from that at the bound- 
ary then a heat flow is produced from the boundary to 
the bulk. Since the bulk is a big reservoir then in the 
main approximation its temperature can be treated as 
time- independent. In this case statistics of becomes 
quasistationary at the second stage (all the correlation 
functions are independent of time). The equations for 
the passive scalar correlation functions are determined 
by Eq. (jl.9JI where the time derivative can be omitted. 
Particularly, at q ^> ru one obtains the closed equation 
for the moments (|2.12l) (where time derivative has to be 
omitted) leading to the same law (|2.19|l . Inside the dif- 
fusion boundary layer ~ where i?o is the tempera- 
ture at the boundary (remind that we assume — in 



the bulk). Then, as it follows from the equation l|l.ll(l . 
(6 n ) ~ ^{m/qf. 

We established that if q ^> 6 at the first stage or 
if q ^> ru at the second stage then the law H2.19JI is 
valid. Moreover, it is possible to formulate the estimates 
(6>") - (9) n (q/5)^ n - 1 > for the first stage and the inequal- 
ities (0 n ) > (0) n (q/ru) 3in ~ 1} for the second stage. The 
estimates show that the high moments of the passive 
scalar are much larger than their Gaussian evaluation, 
this property implies strong intcrmittency in the system. 
The expression (|2.19|l can be treated as a manifestation 
of anomalous scaling. It is extreme since the exponent 
characterizing q-dependence of the moments is indepen- 
dent of n. 



C. Pair Correlation Function 

Here we examine the pair correlation function 
F( r i, r 2)- One could anticipate its non-trivial scaling 
behavior in the region q\,q 2 3> 5 (for the first stage) 
or qi , q<x ^> ru (for the second stage) . If q\ 3> q 2 then 
the only term Vi a [D af 3(ri, r^VipF] in Eq. Ijl.7j) should 
be kept, and we come to the conclusion that F cx q^ 3 . 
Therefore we obtain the behavior, which is similar to the 
passive scalar moments in this asymptotic case. 

Now we consider the case of close points ri and r 2 , 
assuming small value of g^ = ri ll — r 2ll , admitting 
the expressions l|2.3[) and q <§C Q, where Q = (qi + 92) /2 
and q = q\ — qi- Then one obtains from Eq. ljl.7|) 

d t F = fid Q {Q 4 d Q F) + (k/2)8qF + 2nd 2 F 

+ {^Q 3 qd Q d q - 2 M Q 3 gd Q d e 

+ [H w q 2 /g - 2(2 M + 3H 20 )Q 2 g] 8 8 }f 

+ {^iQ 2 q 2 d 2 q - 2H 30 Q 4 g 2 d 2 - 2 f iQ 2 gqd q d e 

+ (2 M + H 20 )Q 2 g 2 d 2 g + H 10 q 2 d 2 e }F (2.20) 

Details of the derivation of Eq. I|2.20|l can be found in 
Appendix |X] It is natural to write F = {0 2 (Q)) (1 + ?) 
Here unity represents the main contribution to the pair 
correlation function, related to the second-order moment, 
and ? is a small correction tending to zero at r± — > r%. 

In the region Q S for the first stage or in the region 
Q ^> ru for the second stage the time derivative in Eq. 
(|2.20() can be omitted in comparison with the turbulent 
diffusion term, which is the first one in the right-hand 
side of Eq. H2.2()|l . If the molecular diffusion terms in 
Eq. (|2.20|) are also neglected then one finds 

{(2H 30 g 2 - 7x 2 )dl - llxd x + 3 + 

[(2 + H 20 )g 2 + H 10 x 2 ]d 2 g + 

(2(1 - 3H 2Q )g + H w x 2 /g) = 0, (2.21) 
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where x — q/Q and we have taken into account that 
{9 2 {Q)) oc Q~ 3 . A solution of the equation l|2.21[l can be 
written as a sum of q, = (qQ) b fb(x, g) where ft satisfy 
the equation H2.21fl (since qQ is zero mode of the operator 
figuring in this equation). It is clear that the main con- 
tribution is related to b = (since negative b are forbid- 
den) . Since the operator in Eq. I|2.21[l has definite scaling 
properties, solutions of the equation can be written in a 
simple self-similar form <; = {q/Q) a ^{Qg/ q)- The prin- 
cipal contribution to c~ is associated with the smallest a. 
Unfortunately, the exponent a is non-universal being de- 
pendent on the coefficients Hy. The molecular diffusion 
smoothes the function ? at q ~ r^/Q. 

The above expressions demonstrate the following nat- 
ural behavior of the pair correlation function. The main 
anomalous dependence of the function is related to the 
second moment of the passive scalar. As to the depen- 
dence on the relative separation between the points r± 
and r2, it is described in terms of the function ? pos- 
sessing a non-trivial scaling behavior with a self-similar 
factor depending on the combination Qg/q. 

III. DECAY ALONG PIPE 

Here we discuss the case when the passive scalar de- 
cays along a pipe in a statistically homogeneous flow, 
which is assumed to be chaotic and has an average veloc- 
ity u along the pipe. Such setup was used by Groisman 
and Steinberg in their experiments fl4j . In this case the 
chaotic flow (elastic turbulence) was excited in a polymer 
solution pushed through a curvilinear pipe. The scalar 
dynamics is then governed by the equation 

d t 9 + ud z 9 + v a \7 a 9 = kV 2 9, (3.1) 

where v is fluctuating part of the velocity (with zero 
mean) and z is coordinate along the pipe. Below, the 
average velocity u is assumed to be much larger, than v 
(that corresponds to the experimental situation). 

If the pressure difference, pushing the flow, is constant, 
the flow is statistically stationary and homogeneous along 
the pipe (which is assumed to have a constant cross- 
section). Therefore u is independent of z and the velocity 
correlation functions do depend on the time differences 
and coordinate differences along the pipe. Thus, for a 
problem with a stationary scalar injection to the pipe 
the scalar statistics is time-independent and the coordi- 
nate z along the pipe plays the role of time in the decay 
problem. Following the procedure described in the first 
section one obtains from Eq. i|3.1|) the following equa- 
tions for the scalar correlation functions 

n n 
ud z F n = K ^ V m^n+ ^2 Vma [ D OLp{ r rn,r k )V kf 3] F n . 
m— 1 m,k— 1 

. (3 " 2) 

The only difference in comparison with Eq. I|1.9|l is that 
the time derivative is substituted by the advection term 
along the pipe. 



As previously, we introduce the coordinate q measur- 
ing a separation of a given point from the wall. The 
average velocity u is a function of q, tending to zero as 
u oc q near the boundary. That leads to the main dif- 
ference between this situation and the one discussed in 
the previous section. The equation for the average scalar 
turns to 

s qd z (9) = [nd^d, + nd 2 ] (9) . (3.3) 

Despite the additional factor q in the left-hand side of 
the equation in comparison with Eq. I|2.5I) . the quali- 
tative picture of scalar evolution remains the same. At 
the first stage the scalar is mostly situated in the layer 
of the width 5 = sq/{ijlz), and the molecular diffusion 
can be neglected. The scalar decay at this stage is alge- 
braic with the longitudinal coordinate z. When 8 reaches 
the boundary layer width ru , the molecular diffusion be- 
comes relevant, and the scalar decay starts to be expo- 
nential. 

As in the previous case, it is possible to obtain com- 
plete statistical properties of the scalar at the first stage. 
In this case one can omit the molecular diffusion term. 
Hence, averages of any single-point functions $(0), such 
as the moments (6 n ) or the scalar PDF P(9,q,z) = 
(S[8 — 6{q 1 z)\) are described by the same equation 

d z ($) - -^ 9 4 cy$) = (3.4) 
soq 

This is a linear equation, and it can be solved using the 
Green function formalism. Namely, the solution of Eq. 
(I3.4|l can be written as 

roc 

(*(*,«))=/ dq'G(z,q,q'Mz = 0,q'), (3.5) 
Jo 

where G is the Green function. In order to obtain an 
explicit equation for the Green function one should first 
solve the corresponding eigenvalue problems. Since the 
operator H is not Hermitian, then one should find both 
it's right and left eigen functions. They are solutions 
of the equations Hf\ + \ f\ = and H + g\ + \g\ = 0, 
respectively, where H + = (so/ ^dqq^dqq^ 1 . We find 

fx = x 3 J 3 (x), g\ = xJ 3 (x), x = 2y/\Jq (3.6) 
Using the orthogonality relation 

/•oc 

/ dqg x (q)f l/ (q) = l6\ 2 5(\-v), 
Jo 

one derives 

/>oo 1 \ 

G(z,q,q') = J o ^exp(-XS- 1 )fx(q)gx(q') 

= q V2 5 {q ,y/2 ex P - <W) h (2-5/ v^ 7 ) • (3.7) 

Since we assumed, that 6 — 9 (q) initially, then for the 
n-th order scalar moment {9 n {q,z)) the initial condition 
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is 0£(q). If 6 <C L then one can substitute (0 n (Q, z)) by 
where $o = #o(0,0). In this case an explicit integra- 
tion in Eq. (|3.5() leads to the expression 

(9 n {z, q)) = ^ ex P & (1, 4, J/g) . (3.8) 

Again, we obtain a universal profile. If q 3> 8 then 
both the exponent and \F\ can be substituted by unity 
to obtain (6 n (z,q)) = ^S 3 /(6q 3 ). If q < 6 then 
(9 n (z,q)) = i?q . In order to obtain the scalar PDF, 
one should use the initial distribution P(z = 0, q, 9) = 
6(9 - o (g)) = -(l/^(g))% - q (9)). Here g o (0) is de- 
fined, as previously, by the relation #0(90) = The 
convolution with the Green function is now 



P(z,q,9) = / dq'G(z,q,q')P(0,q,9) 



(3.9) 



q 3/2 q 



1/2 



■ exp 



q qo 



At large distances, g > <5 2 /i (and also L ^> 5) one can 
use the approximation: 



6g 3 g o 2 o (g o ) 



exp (-5/g) 



(3.10) 



Again this expression cannot be exploited for calculation 
of moments of the passive scalar because of the diver- 
gency near the maximum value of 9, $0. 

For the pipe, the characteristic length where the dif- 
fusion becomes relevant is determined by the same ex- 
pression 12.41) . When 8 diminishes down to ru another 
regime comes. For the density of pollutants the regime 
is characterized by an exponential decay of the passive 
scalar moments along the pipe. For the first moment the 
decrement of the decay can be extracted from Eq. 13.311 , 
it is determined by the smallest eigen value of the opera- 
tor in the right-hand side of the equation. Thus, for large 
enough z the decay law (9) oc exp(— az) is correct where 
a ~ K 1 / 4 /i 3 / 4 s " 1 . The eigenvalue problem, correspond- 
ing to the equation (|3.3|l can be solved numerically, then 
one obtains a « 3.72k 1 / 4 /x 3 / 4 s " 1 . 

For the temperature, fixed at the boundary, we ob- 
tain a quasi-stationary distribution if the fixed value is 
z-independent. The average temperature is described by 
the equation I3.3fl where z-derivative is dropped. Then 
we return to the same equation as previously (for the 
decay in the vessel), with the solution l|2.11|l . Higher 
moments of the temperature are slightly different than 
for the decay in the vessel. However, qualitatively their 
behavior is the same. Namely, (9 n ) w $ l if 1 ^ r bi, an d 
(6 n ) ~i?5(r M /?) 3 , if g»r &i . 

One can examine the pair correlation function F of 
the passive scalar. Qualitatively, it is the same, as for 
the decay in the vessel. Namely, the main anomalous de- 
pendence of the function is related to the second moment 
of the passive scalar. The function is practically undis- 
tinguished from the second moment if gi , g2 <C r&j . As 



to the dependence on the relative separation between the 
points 7*1 and r-i for gi, qi 3> ru, it is described in terms 
of a function Q, (# 2 (Q)) (1 + possessing a non-trivial 
scaling behavior with a self-similar factor depending on 
the combination Qg/q. 



IV. CONCLUSION 

We have investigated the passive scalar (concentration 
of pollutants or temperature) evolution in the chaotic and 
turbulent flows near the boundary (wall), which domi- 
nates the advanced stages of the passive scalar homoge- 
nization. There are some universal features of the decay 
related to the universal velocity dependence on the coor- 
dinate perpendicular to the wall. We assumed that the 
Peclet number is large enough to guarantee the condition 
L 3> ru where L is the width of the peripheral region 
and ru is the thickness of the boundary layer (which is 
determined by equating the molecular diffusion and the 
turbulent diffusion). The width of the peripheral region 
is the thickness of the viscous boundary layer for the de- 
veloped high-Reynolds turbulence and is of the order of 
the vessel size for the elastic turbulence. For the decay 
along a pipe L is of the order of the pipe radius. 

We have considered the evolution beginning after that 
homogenization of the passive scalar is finished in the 
bulk. Then the initial passive scalar distribution is char- 
acterized by the length L, which is much larger than ru- 
The first stage of the passive scalar decay in the periph- 
eral region is insensitive to diffusion. The characteristic 
length of the passive scalar distribution at this stage <5 
satisfies the inequalities L ^> 8 ru- For the decay in 
a closed vessel the length behaves as 8 oc t~ x / 2 , whereas 
for the decay along pipe 8 cx z _1 (where z is coordi- 
nate along the pipe). The passive scalar inhomogene- 
ity is kept mainly for separations g from the boundary, 
satisfying q < 8 whereas for larger q the advection effec- 
tively homogenize the passive scalar (the mechanism is in 
stretching producing small-scale fluctuations of the pas- 
sive scalar, which are effectively removed by diffusion). 

When 8 achieves the thickness of the diffusion bound- 
ary layer, ru, the decay of the concentration of pollutants 
becomes exponential. This character of the decay is ex- 
plained by the space distribution of the passive scalar 
inhomogeneities, which are concentrated mainly in the 
diffusive boundary layer, and decay due to the flow to 
the bulk, which is proportional to the level of the inhomo- 
geneities. The decrements of the decay can be expressed 
in terms of the Peclet number which is Pe = VlL/k. 
For the viscous boundary layer the Peclet number coin- 
cides with the Schmidt number, whereas for the elastic 
turbulence it differs from the Schmidt number by the fac- 
tor which is the Reynolds number. For the time decay 
the decrement is proportional to Pe ' , whereas for the 
decay along the pipe the decrement is proportional to 
Pe" 1 / 4 . For the temperature fixed at the boundary the 
second stage is characterized by a quasi-stationary dis- 
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tribution, since the bulk serves as a big reservoir keeping 
its temperature constant and absorbing heat, transferred 
from the boundary. 

Our theoretical predictions can be compared with the 
experimental data [lj, Il9f . obtained for the dye con- 
centration decay (homogenization) along the curvilinear 
pipe, where the chaotic flow (elastic turbulence) is excited 
in a dilute polymer solution. The most interesting com- 
parison concerns the Pe-dependence of different quanti- 
ties, which was extracted from the experiment |19| by 
varying the diffusion coefficient k by two orders of mag- 
nitude (it was done by changing the molecular weight of 
the dye carriers). The experimental data (correspond- 
ing to the second stage in our theoretical scheme) are 
in a good agreement with the law Pe -1 ^ 4 both for the 
thickness of the diffusion boundary layer, ru, (which is 
determined experimentally by a maximum in the passive 
scalar fluctuations) and for the decrements of the passive 
scalar decay along the pipe. 

The peripheral region serves as a source supplying the 
passive scalar to the bulk. Since the passive scalar decay 
in the peripheral region is slow (comparing to the bulk), 
it can be considered as a quasi-stationary source of the 
passive scalar for the bulk, where the passive scalar cor- 
relation functions adjust adiabatically to the level of the 
supply. Therefore the situation is close to one charac- 
teristic of the passive scalar, supplied by pumping, in a 
smooth random velocity field. Statistical properties of 
the passive scalar in this case were examined in the pa- 
pers pf |8|, |9( . Particularly, correlation functions of the 
passive scalar in the bulk should possess a behavior close 



to the logarithmic one. This prediction is also in agree- 
ment with the experimental data 0, 0] . 

One of remarkable predictions of our theory is strong 
intermittency of the passive scalar and extreme anoma- 
lous scaling of its moments for the dependence on the 
separation from the boundary q: all the moments scale 
identically as q~ 3 for L ^> q ^> 8 or L ^> q ^> ru, de- 
pending on the stage. Physically that means that the 
passive scalar decay is related to rare events like escapes 
of separate passive scalar blobs from the boundary layer 
with their subsequent transportation to the bulk. As to 
the correlation functions of the passive scalar, our theory 
predicts some definite self-similar and scaling behavior, 
which are complex due to the space inhomogeneity char- 
acteristic of the problem. It would be interesting to check 
the predictions experimentally. 

The final remark concerns an extension of our results 
to other problems. As it was noted in the paper Q, 
the scheme developed for the passive scalar decay can 
be without serious modifications applied to fast binary 
chemical reactions. We believe that minor modifications 
of the scheme can make it applicable for more compli- 
cated chemical reactions. The other problem, which can 
be posed for the peripheral region, is dynamics of poly- 
mers (say, for the case of the elastic turbulence). Then 
one should go away the scope of the passive approach. 
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APPENDIX A: PAIR CORRELATION FUNCTION 



The anomalous diffusion operator A in the dynamics of pair correlation functions given by the last terms of the 
equation i|1.7fl . It can be rewritten as 



A = D l - 



d d 



(3D 1 ' 



d 



V dX l / dXJ 



dX l dXi 

Here X % is extended set of coordinates X = {q\, ru, ri2, q 2 , T21, ^22}, and 

H 3 q 2 q 2 



(Al) 



D 



Ml 






H 3 qh 2 2 



\ liq\q2{r 12 - ^22) 



H w q\ M\q\{ r 2\ - ru) 

H la q\ iiqiql(r 2 2 - r 12 ) 

mq 2 { r 2 2 - r 12 ) Liqiql(r 2 i - r n ) 



Hiqiq 2 + d 1 
d 21 



H x qxq 2 +d 2 
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Mi92(n.i - r 2 i) liq\q 2 {r 12 - r 22 ) \ 

19192 + d 11 d 12 

d 21 H lQl q 2 + d 22 





Hwq 2 






Hioql 



H, = H 
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(At + 3H 20 /2) ((r u - r 2 i) 2 + (r 12 - r 22 ) 2 ) 
H 3 =n + H 30 ((r u - r 2 i) 2 + (r 12 - r 22 ) 2 ) 
d l] = H 20 qiq 2 (ru - r-n){r\j - r 2 j) 
Now we turn to the new variable set Y = {Q,q,Ri,R 2 , g,(j)}, such that Q — {q% +q 2 )/2,q 



(A2) 
(A3) 
(A4) 
(A5) 

= q%- q 2 ,Ri = (ru + 



r2i)/2, gcoscf) = ru — r2i, psm</> = r\ 2 — r 22 . In this case, the anomalous diffusion operator transforms to 



,1 /)'"//• 



_d d_ 

3 dY k gyl 



JJIJ rpk 





d 




\dY k j 


dY l + 


\ dX* ) 



rp k 







Q Y k 



_d d_ 



iJL 
dY* 



(A6) 
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where Tj = dY l /QX 7 is the Jacoby matrix: 



T = 



/1/2 








1/2 





° \ 


1 








-1 











1/2 








1/2 











1/2 








1/2 





COS</> 


sin </> 





— COS0 


— sin (f> 


V o 


-sin((/»)/e 


cos((j))/g 





shiOVe 


- COs((j))/Q J 



(A7) 



After some algebraic manipulations one can easily obtain the new diffusion tensor and the first order term /*. 
The exact expressions are rather bulky. However, if the initial distribution of scalar depends only on q, the only 
relevant terms, are those, preceding dQ,d qi d e . Therefore one can reduce the final variables set to {Q,q, g}- In the 
limit Q » q, p — > one has 



Q\p + H 30 g 2 /2) 
D= \ 2pQ 3 q 
-pQ 3 g 



2pQ 3 q -^Q 3 e 

4/iq 2 Q 2 - 2H 30 Q i g 2 -pQ 2 qg 

-nQ 2 qg (2p + H 2a )Q 2 + H w q 2 



(A8) 



/« = 4Q 3 (p + H 30 p 2 /2), P = 0, 



fP= HW__ 4{fi + 3H2o/2)Q 2 p 



(A9) 



This yields the following expression for A: 



A=(n + H 30 p 2 /2)d Q Q i d Q + 4nQ s qd Q d q - 2pQ 3 gd Q d e + 
+ (4pQ 2 q 2 - 2H 3Q Q 4 g 2 )d 2 - 2pQ 2 qgd q d e + [{2p + H 2Q )Q 2 g 2 + H 1Q q 2 ] d 2 g + 



+4Q 3 (p + H 30 p 2 /2)d Q + 



H w q 2 



4( M + 3H 20 /2)Q 2 p 



(A10) 
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